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ABSTRACT 

The large range of time and length scales involved in type la supernovae (SN 
la) requires the use of flame models. As a prelude to exploring various options 
for flame models, we consider, in this paper, high-resolution three-dimensional 
simulations of the small-scale dynamics of nuclear flames in the supernova envi- 
ronment in which the details of the flame structure are fully resolved. The range 
of densities examined, 1 to 8 x 10 7 g cm -3 , spans the transition from the laminar 
flamelet regime to the distributed burning regime where small scale turbulence 
disrupts the flame. The use of a low Mach number algorithm facilitates the ac- 
curate resolution of the thermal structure of the flame and the inviscid turbulent 
kinetic energy cascade, while implicitly incorporating kinetic energy dissipation 
at the grid-scale cutoff. For an assumed background of isotropic Kolmogorov 
turbulence with an energy characteristic of SN la, we find a transition density 
between 1 and 3 x 10 7 g cm -3 where the nature of the burning changes qual- 
itatively. By 1 x 10 7 g cm -3 , energy diffusion by conduction and radiation is 
exceeded, on the flame scale, by turbulent advection. As a result, the effective 
Lewis Number approaches unity. That is, the flame resembles a laminar flame, 
but is turbulently broadened with an effective diffusion coefficient, ~ u'l, 
where u' is the turbulent intensity and / is the integral scale. For the larger 
integral scales characteristic of a real supernova, the flame structure is predicted 
to become complex and unsteady. Implications for a possible transition to deto- 
nation are discussed. 

Subject headings: supernovae: general — white dwarfs — hydrodynamics - 
nuclear reactions, nucleosynthesis, abundances — conduction — methods: nu- 
merical — turbulence 



Lawrence Berkeley National Laboratory, 1 Cyclotron Road, MS 50A-1148, Berkeley, CA 94720 

2 Dcpartment of Astronomy and Astrophysics, University of California at Santa Cruz, Santa Cruz, CA 
95064 

3 Dcpartment of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794 



-2- 



INTRODUCTION 



The complex small-scale dynamics of turbulent thermonuclear flames are essential to 
understanding SN la. Since the range of length scales is so large (star size ~ 10 8 cm; 
Kolmogorov scale and flame thickness < 1 mm), full star calculations must use a sub- 
grid model to describe t he b urning on the unresolved scales. Turbulent flame models 
( INiemeyer fc Hillebrandtl 119951 1 th at move the flame at the turbulent velocity have been 
successful at producing explo sions (iRopke fc Hillebrandtl 120051 1 with pure deflagrations. An- 
other popular flame model (jKhokhlovl Il99ll ) moves the flame at the spee d dictated by 



the R ayleigh- Taylor instability, and has also produced successful explosions (IGamezo et al. 
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20031 ). although the authors argue tha t a more realistic sup ernova is produced when a det- 
onation ensues late in the explosion ( Gamezo et al.l 120051 ) . These two flame models can 



produce large differences in the flame speeds at small scales, which will have a large effect 
on the outcome of full-star simulations. Validation of these models against resolved flame 
calculations is needed to resolve this discrepancy. 

As the flame propagates outward from the center of the star, it encounters lower den- 
sities, the speed decreases, and the flame thickens. A critical length scale in turbulent 
combustion is the Gibson scale - the l ength scale a t which the laminar flame speed equals 
the turbulent velocity (see for example |PetersH2000l ). 
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Here sl is the laminar flame speed, and u 1 and / are velocity and integral length scales 
characterizing the turbulence, and we have assumed Kolmogorov scaling. At a density of 
around 10 7 g cm -3 , the flame be comes thick enough that t urbulent eddies can disrupt its 



structure before they burn away (INiemeyer fc Woosleylll997l ). that is, the flame thickness is 
larger than the Gibson scale. At this point, the burning fundament ally changes character 
and the flame is said to be in the distributed burning regime (see iPetersI |2000| for some 
discussion) . 

Here we look at the interaction of the flame and turbulence on the scale of the flame, 
with the aim of validating turbulent flame models and better understanding the transi- 
tion to distributed burning. Previous studies of flames interacting with turbulence have 
focused on flames interacting with vortical flow (two-dimensional simulatio ns presented in 



2000a 



Ropke et al.l 120041 ) or statistical methods using one-dimensional turbulence (ILisewski et al. 



b). An open question is whether in the distributed burning regime, a mixed region 



of partially burned fuel and ash can grow large enough such that it can ignite a detonation 
( INiemeyer fc Woosleylll997l ; iKhokhlov et al.lll997l ; lNiemeyerlll999l ) . Detonations have gained 
renewed interest lately as a means to burn the carbon/oxygen fuel left behind near the center 
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of the white dwarf in pure deflagration models (IGamezo et al.ll2005l ). Although, it is unclear 
whet her detonations can traverse the pockets of partially burned fuel (IGolombek fc Niemeyer 
20051 ). Recent results also suggest that the size of the reg ion of partially burned f uel needed 



to init iate a detonation is larger than previously believed fjDursi fc Timmesll2006l ). iPan et al. 



( 120081 ) discusses the role of turbulent intermittency on the conditions needed for transition 
to detonation. 

Previously, w e showed the transition to the distributed burning regime in two dimensions 



(IBell et al.l l2004bl ). where the Rayleigh- Taylor instability (RT), growing on scales smaller 
than the flame thickness it self, is responsible for the transition. It has been suggested 



(INiemeyer fc Kersteinlll997l ) that turbulence generated by RT in type la supernovae should 



obey Bolgiano-Obukhov (BO) statistics - formulated from considering a potential energy 



' Chertkov 


2003; 


Zineale et al. 


2005a) 



have shown that turbulence is indeed Kolmogorov in nature and becomes isotropic on the 
small scales ( iZingale et al. I l2005bl : ICabot &: CoShood ) . As first discussed in lNiemeyer fc Kersteinl 



(119971 ). BO scaling leads to a lower transition to distributed burning, and therefore, makes 



a detonation transition more difficult. The refore, although o ur two-dimensional results led 



us to conclude that a DDT was unlikely (IBell et al.l l2004bl ). a three-dimensional study is 
warranted. 

Figure [TJ shows how the relevant length scales vary with fuel density. In particular, the 
two red curves compare the laminar flame width and the Gibson scale. The laminar flame 
widths have been calculated from simulated one-dimensional profiles. The Gibson scale is 
ev aluated assuming a tu r bulen t intensity, u* , of 10 7 cm s" 1 on a length scale, L*, of 10 6 cm, as 



in 



Niemeyer fc Woosleyl (119971 ). We see that the transition to distributed burning, assuming 
Kolmogorov turbulence, should occur around 3 x 10 7 g cm -3 . 

In this paper, we present a three-dimensional study of turbulent thermonuclear flames 
that explores a range of conditions from the flamelet regime to the distributed burning 
regime. Specifically, we use a flame sheet embedded in a maintained turbulent velocity field. 
Of particular interest is the response of the flames to wrinkling by the turbulence, and the 
scaling of the turbulent flame speed with the turbulent intensity. Therm al diffusion is orders 
of ma gnitude greater than species diffusion (the Lewis number is large) (ITimmes fc Woosley 
19921 ) and has a significant effect on the turbulence-flame interactions. Attention is then 
focused on the distributed burning regime and potential implications for a transition to 
detonation are discussed. 
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SIMULATION DESCRIPTION 



We use a low Mac h number hydrody namics code, adapted to the study of thermonuclear 
flames, as described in lBell et al.l (j2004al ). The advantage of this method is that sound waves 
are filtered out analytically, so the time step is set by the the bulk fluid velocity and not the 
sound speed. This is an enormous efficiency gain for low speed flames. The input physics 
used in the present simulations is largely uncha nged, with the exce ption of the addition of 
Coulomb screening;, taken from the Kepler code flWeaver et al.lll978f ) to the 12 C( 12 C,7) 24 Mg 
reaction rate. This yields a small enhancement to t he flame speed , and is included for 
completeness. The conductivities are those reported in iTimmesi ( 120001 ) . and the equation of 
state is the Helmholtz free-energy based general stellar EOS described in iTimmes &: Swesty 
( 120001 ). We note that we do not utilize the Coulomb corrections to the electron gas in the 
general EOS, as these are expected to be minor in the conditions considered. 

The non-oscillatory finite-volume scheme employed here permits the use of implicit 
large eddy simulation (iles). This technique captures the inviscid cascade of kinetic energy 
through the inertial range, while the numerical error acts in a way that emulates the physical 
effects of the dynamics at the grid scale, without t he expense of reso lving the entire dissi- 
pation subrange. The appr oach was introduced by Boris et al. (fl992), and has sin c e bee n 
use d by many authors (e.g 
andlMar golin et aLJ 2006|) 



Youngj (ll99lh . |Porter et al.l (119921 ) . iGrinstein & Furebv! (l2004h. 



An overview of the technique can be found in IGrinstein et al. 



( 120071 ). I Asp den et al.l (120081 ) presented a detailed study of the technique using the present 



numerical scheme, including a characterization that allowed for an effective viscosity to be 
derived. Thermal diffusion plays a significant role in the flame dynamics, and so is explic- 
itly included in the model, whereas species diffusion is significantly smaller, and so is not 
explicitly included. 

The tur bulent velocity field w as maintained using a forcing term similar to that used in 
the study of I Asp den et al.l (120081 ) . Specifically, a forcing term was included in the momen- 
tum equations consisting of a superposition of long wavelength Fourier modes with random 
amplitudes and phases. The forcing term is scaled by density so that the forcing is somewhat 
reduced in the ash. This approach provides a way to embed the flame in a zero-mean turbu- 
lent background, mimicking the much larger inertial range that these flames would experience 
in a type la supernova, without the need to resolve the large-scale convective motions that 
drive the turbulent e nergy cascade. Figur e [2] shows an example kinetic energy wavenumber 
spectrum taken from I Asp den et al.l (120081 ) at a resolution comparable to that used here; the 
dashed line denote s a minus five - thirds decay, and is illustrative of the turbulence found in 
the present study. lAspden et al.l (120081 ) demonstrated that the effective Kolmogorov length 
scale is approximately 0.28Ax, and the integral length scale is approximately a tenth of the 
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domain width. 

Figure [3] shows the simulation setup. The simulations were initialized with carbon fuel 
in the lower part of the domain and magnesium ash in the upper, resulting in a downward 
propagating flame. A high-aspect ratio domain was used to allow the flame sufficient space 
to propagate. Periodic boundary conditions were prescribed laterally, along with a free-slip 
base, and outflow at the upper boundary. 

Five simulations were run to investigate the turbulence-flame interactions in different 
burning regimes. The five cases will be referred to as cases (a) through (e). The study was 
designed so that Ig/Il varied across a number of orders of magnitude, from approximately 
4 x 10 3 in case (a) corresponding to the flamelet regime, to 2 x 10 -5 in case (e) corresponding 
to the distributed burning regime. The different burning regimes were achieved by varying 
the density of the carbon fuel from 8 x 10 7 for case (a) down to f x fO 7 in case (e). Case 
(a) corresponds to near the center of a supernova, has a thin laminar flame and relatively- 
low turbulence levels. Case (e) corresponds to the conditions near the edge of a supernova, 
has a thicker laminar flame and relatively more intense turbulence. The derivation of the 
turbulence conditions will be described below. 

For each case, a flat laminar flame simulation was first run to steady-state to establish a 
laminar flame width and speed. The laminar flame width was taken to be the carbon width, 
specifically II = (AXc)/ max|VXc|, where Xq is the carbon mass fraction. The variation 
of laminar flame thickness with fuel density is shown by the solid red curve in figure [TJ The 
laminar flame solution also provides initial flame conditions for the turbulent simulations. 
Figure H] shows normalized laminar profiles for the five cases. The length scale has been 
normalized by the carbon width, and each quantity has been normalized by the minimum 
and maximum value across the domain, i.e. q = (q — 5 m in)/(?max _ ?min) for a generic quantity 
q. Properties of the laminar flames are summarized in Table [TJ 

Based on the resolution studies presented in iBell et al. (bo04al lbh. the resolution for the 



turbulent flame simulations was chosen to have four cells across the laminar flame width 
Specifically, the carbon profile from each laminar flame simulation was mea sured, and 



Ax se t to Zl/4, which is shown by the solid blue line in figure [U and, following I Asp den et al. 



(120081 ). the effective Kolmogorov length scale is shown by the solid green line. This choice 
of cell width corresponds to many more cells across the entire flame; at the resolution of the 
turbulent flame simulations to follow, the spatial extent shown in figure H] would be spanned 
by twenty computational cells. Since it is expected that the turbulence will thicken the 
flame, this choice of resolution is believed to be sufficient to resolve the reaction zone of the 
turbulent flame. 
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Computational expense restricts the domain to a cross-section of 256 cells, and so, since 
Ax is known, the domain width follows, and the forcing term determines the integral length 
scale. The domain width and integral length scales are shown in figure [TJ by the solid and 
dashed black lines, respectively. It should be emphasized that the integral length scale 
and effective Kolmogorov length scales are restricted by computational expense and are not 
reflective of the values found in a supernova. 

Turbulence characteristics can be estimated using a simple representation shear gener- 
ation by the Rayleigh- Taylor instability. On the large scale of the supernova, the flame will 
be unstable with a typical feature size of L* ~ 10 6 cm . The velocity on this scale can be 



estimated using the Sharp- Wheeler model (jSharplll984l ) 



u 



(2) 

For an Atwood number, At ~ 0.5 and gravitational accelerat ion, q ~ 10 9 cm s~ 2 , one has 



u*(L*) ~ 10 7 cm s 1 , in agre ement with t he numbers used in iNiemeyer fc Woosleyl (119971 ). 



More detailed simulations by iRopkd (120071 ) in 3D find a distribution of turbulent intensities 
peaking at 10 7 cm s _1 , but with a high-velocity tail extending out to 10 8 cm s -1 . For the 
current simulations, the turbulent kinetic energy released by the large-scale RT motions is 
assumed to cascade down to the flame scale following a Kolmogorov spectrum, as 



u\l) 



«W ( 77 



1/3 



(3) 



therefore, given the integral length scale, I, the turbulent velocity, u'(l), can be computed. 
Throughout this paper, u* = 10 7 cm s _1 and L* = 10 6 cm have been assumed. 

In cases (a)-(d), the aspect ratio of the domain was 4:1 and the flame placed one do- 
main width from the upper boundary. In case (e) an aspect ratio of 8:1 was required to 
accommodate the flame in the distributed burning regime, and the flame was initially placed 
at the midpoint of the domain. In each case, adaptive mesh refinement was used to reduce 
the computational expense. A base grid with 128 cells across was used, with one level of 
AMR with refinement factor 2 concentrated around the flame sheet. This gave an effective 
resolution of 256 x 256 x 1024 for cases (a)-(d) and 256 x 256 x 2048 for case (e). Experi- 
ments with additional resolution confirmed the adequacy of this choice. It was also ensured 
that the use of AMR had no detrimental effect on the tur bulence. Spec i fically , the energy 
containing scales were well captured by the base grid (see I Asp den et al.l (120081 ) for a study 
of maintained homogeneous isotropic turbulence using the same approach) and there was 
adequate refinement ahead of the flame to ensure that the cascade of energy to smaller scales 
in the refined region (which have a short turn over time) could occur before these scales could 
interact with the flame. Table [2] summarizes the properties of the turbulence simulations. 
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3. RESULTS 

In each case, the flow undergoes an initial transition, as the turbulence wrinkles and 
thickens the flame, until a quasi-steady state is established. Figure [5] shows instantaneous 
vertical slices of fuel consumption rate (left panel) and temperature (right panel) for each 
of the five cases once the quasi-steady state has been reached. In each case, the values have 
been normalized by the corresponding laminar values, with the exception of case (e) where 
the local fuel consumption rate is much lower and so has been normalized by a fifth of the 
laminar value for contrast. By construction, the length scales are normalized by the laminar 
flame width. Half of the domain is shown for cases (a)-(c) and case (e), and the entire 
domain is shown for case (d). 

Case (a) presents smooth and even burning, and is perturbed very little by the back- 
ground turbulence. The temperature profile remains sharp. In case (b), the flame surface has 
been deformed by the turbulence. Both regions of enhanced burning and regions of decreased 
burning are observed, and appear to be correlated with the curvature of the flame sheet. 
Specifically, enhanced burning appears to occur where the center of curvature is within the 
fuel, and decreased when the center of curvature is in the ash. The temperature field presents 
regions that are sharp, and regions that appear to be more diffuse. Again, this appears to 
be correlated with curvature, with the more diffusive regions occurring where the center of 
curvature is in the products. 

In cases (c) and (d), as Ig/Il decreases further, the background turbulence becomes 
increasingly influential; the temperature field becomes more mixed, and the deformation of 
the flame surface increases. The burning appears to occur in small high-intensity pockets, 
punctuated by regions of local extinction. 

In case (e), a dramatically different burning mode is observed. The temperature mixed 
region and the burning region are much broader. The burning appears to be much less 
intense (recall the image has been normalized by a fifth of the laminar value for contrast) 
and is restricted to the high temperature end of the mixing zone. There is no well-defined 
flame surface, but a broad flame brush. Interestingly, there appears to be some residual 
low-level burning well above the main burning region, suggestive of incomplete burning in 
the main flame zone. 

These observations are further reinforced by three-dimensional renderings of the fuel 
consumption rate, shown in figure [6], which elucidates the three-dimensional structure of the 
flames. The images have been scaled in the same way as the corresponding slices. Case (a) 
burns as a coherent connected flame sheet, but as the relative turbulence level increases in 
cases (b) to (d), the flame sheet becomes increasingly disrupted and presents regions of local 
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extinction. Finally in case (e), in the distributed burning regime, the behavior is completely 
different with a flame brush that is much broader but burns less intensely than the laminar 
flame. 

The curvature effects observed and the resulting burning rates are a consequence of the 
thermo-diffusively-stable nature of the flames; the Lewis number is high - thermal diffusion 
is much greater than species diffusion. Where the background turbulence elevates part of 
the flame surface, resulting in negative curvature, there is a focusing of heat by diffusion. 
Consequently, as the reaction rate is extremely sensitive to temperature, the fuel burns 
quickly, and the flame is flattened. Conversely, where the turbulence pushes the flame 
downwards, creating positive curvature, temperature diffusion leads to a defocusing of heat, 
and the burning rate decreases. Again, this tends to flatten the flame. Curvature effects will 
be explored further below. 

The instantaneous global turbulent flame speed normalized by the laminar flame speeds 
are shown in figure the time scale is normalized by the laminar burning time - the time 
it takes the laminar flame to burn one flame width. Note the different time axis for case 
(e). Case (a) burns almost exactly at the the laminar flame speed. Curiously, cases (b) and 
(c) actually burn more slowly overall than their laminar counterparts, and case (d) is only 
around 30% faster. However, case (e) burns between approximately five and six times the 
laminar flame speed. 

To investigate the effects of curvature, it is useful to be able to define a flame surface. 
There are a number of ways to do this, but it has been found that using isosurfaces of 
temperature or fuel mass fraction is a practical way that avoids the difficulties associated 
with the local extinction in these particular flames. We choose values for the isosurfaces 
based on the temperature or fuel mass fraction values at the peak local fuel consumption 
rate in the laminar flame. The isosurface was then located by usin g a standard "m arching 



cubes" algorithm, and optimized using the "QSlim" algorithm, see iGarlandl (119991 ) . Since 
the temperature field is more diffusive than the fuel, the resulting isosurface was found to 
be much smoother for temperature than the carbon, particularly where the isosurface was 
within the fuel itself. Furthermore, because of the higher levels of turbulence in cases (d) 
and (e), the isosurfaces were found to be too contorted to be useful. 

Since the fuel consumption occurs over a finite width and is not localized to the flame 
surface, a local consumption-based flame speed was evaluated by integrating the fuel con- 
sumption rate through the surface. This involved constructing a set of integral curves along 
the gradient vector of the progress variable through each point on the isosurface, ensuring 
that the curves extended well beyond the region where the fuel consumption rate decreases 
to zero. (We note that for this construction, we orient the normal to point into the ash. 
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With this orientation, the curvature, k is negative when the center of curvature is in the 
fuel.) These integral curves provided bounding edges of prism-shaped subvolumes that ef- 
fectively cover the reaction zone. The consumption-based flame propagation speed was then 
computed over each subvolume of the reaction zone by integrating the computed fuel con- 
sumption rate over the subvolume. The local consumption-based flame speed was defined 
as 

4= y A \ p^Fdn (4) 

where poX F)0 is the initial fuel density, fl is the subvolume, pujp is the local fuel consumption 
rate, and A is the area of intersection of the flame with Q. Defining the local speed in this 
way has the property that the global burning speed is the integ ral of s l T oyer the isosurface. 



For additional detail about construction of the elements, Q, see iBell et al.l ( 120051 ) . 



Figure [8] compares instantaneous temperature and carbon isosurfaces for case (b), and 
is colored by the integrated fuel consumption rate, s l T . The correlation between burning rate 
and curvature observed in figure [5] is more apparent here. Specifically, the regions of high fuel 
consumption rate are strongly correlated with regions of negative curvature, i.e. where the 
centers of curvature are within the fuel. Two examples are highlighted by arrows labeled 'H'. 
This correlation appears to differ slightly between the two isosurfaces. On the temperature 
surface, there are regions of negative curvature where the burning is very low; two examples 
are shown by arrows labeled 'L'. However, the carbon surface does not appear to pass through 
these regions; the labels are in the same place, but the carbon surface is much lower. The 
temperature and carbon mass fraction have become decorrelated. Regions exists around the 
peak burning temperature where there is no fuel, and so burning cannot occur. This is a 
direct result of the competition between mixing by turbulence and thermal diffusion. 

Figure M shows joint probability density functions (jpdf) of curvature (normalized by 
the laminar flame width) and local integrated fuel consumption rate (normalized by the 
laminar value), based on the temperature isosurface, ensemble-averaged over a number of 
time-points after the flames have reached a quasi-steady state. The (negative) correlation 
is clear and quantifies the between curvature and local flame speed observed above. It is 
also evident that there is a significant change in burning across the three cases. Case (a) 
shows that the majority of the fuel consumption occurs around the laminar rate. Case (b) 
demonstrates that the flame is bimodal, in the sense that there are regions that are burning, 
which are negatively correlated with curvature, but there are also regions of low burning or 
local extinction, which occur with curvature of both signs. This suggests that where the 
flame is burning, it is burning close to the laminar value (with low probability variability), 
but because there are also regions of local extinction, the overall burning rate per unit area 
is lower than the laminar value. Finally, in case (c), there is no preferred rate of burning, a 
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slight negative correlation with curvature, and significant regions of local extinction. 

A common approach used for determining a model flame speed involves assuming that 
the flame is burning locally at the laminar flame speed, and that global burning rate is equal 
to the laminar flame speed times the area of the flame. Figure [10] shows the normalized 
burning rate per unit area for cases (a-c). If the flame- modeling approach described above 
was appropriate for these flames, this measure would be close to unity. This is indeed seen 
to be the case for case (a), as the flame is burning in a very similar way to the laminar flame. 
However, it is too simplistic an approximation for other two cases. The strong response 
to curvature modulates the local burning speed. The values of approximately 0.8 and 0.5 
suggest that using the laminar flame speed times the area of the flame would lead to an 
overestimation of the global burning by factors of approximately 1.25 and 2 respectively. 
The analysis indicates that even in the flamelet regime, a Markstein correction is needed to 
predict the local flame speed accurately. 

The way in which the burning occurs can be analyzed further by considering other joint 
probability density functions. Figure [TT] shows the jpdfs for all five cases of temperature and 
carbon density, and figure [T21 shows fuel consumption rate against temperature. Here the first 
moment has been taken with respect to fuel consumption rate, specifically, for jpdf P(Q,T) 
for fuel consumption rate Q and temperature T, the first moment is defined as QP(Q,T). 
This highlights the burning regions and prevent the large proportion of the domain that 
is not burning from dominating the jpdf. Again, the PDFs have been ensemble-averaged. 
The solid red line in each case denotes the laminar flame correlation. In case (a), it is clear 
that the burning follows an almost identical path to the laminar flame. Case (b) present a 
similar correlation, but there is greater variability around the main burning path. In cases 
(c) and (d), the main burning path appears to have shifted away from the laminar flame, and 
there is a large amount of variability from that main path. In case (e), a dramatic change is 
observed, the burning path has collapsed to a single path that is significantly removed from 
the laminar flame. In this case, the turbulent mixing dominates the thermal diffusion and 
so the temperature and fuel cannot become decorrelated. Therefore, a single burning path 
is observed, which is different than the laminar flame, with very little variability. 

Figure [T3] shows the joint probability density function of fuel consumption rate and 
carbon mass fraction for case (e), again the first moment with respect to fuel consumption 
rate has been taken to highlight the burning region. Similar to the corresponding temperature 
plot, the burning path has collapsed to a curve that is significantly different than the laminar 
flame and there is little variability around it. Importantly, this figure demonstrates that the 
burning occurs at much lower mass fractions than in the laminar flame. 

A normalized instantaneous flame structure in case (e) is compared with the laminar 
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flame structure in figure [T4J The turbulent flame structure has been obtained by taking 
planar averages. Here, the length scale has been normalized by the laminar flame carbon 
width, and each of the other quantities has been normalized by the corresponding laminar 
values, expect for the velocity, which has been normalized by the turbulent burning speed 
due to its significantly enhanced value. The turbulent structure is dramatically different. 
Note how the shape of the temperature and density profiles have changed, becoming closer 
to hyperbolic tangents. In particular, the width of the temperature profile has decreased 
relative to the carbon profile because the mixing is dominated by the turbulence rather 
than thermal diffusion. Therefore, any comparison between the turbulent and laminar flame 
widths will depend on the arbitrary choice of width definition. In spite of this, a rough 
estimate of the carbon width is approximately 60 times the laminar flame width (roughly 140 
cm) . Some effect of thermal diffusion is still evident as the temperature profile and therefore 
density profile are slightly wider than the species profile; the effective Lewis number is still 
greater than unity. One of the biggest differences between the two flame structures is the fuel 
consumption rate. In the turbulent case, the fuel consumption is very much lower, the peak 
value is approximately an eighth of the laminar value. The physical location also appears 
to have shifted to a relatively higher z location. It should be noted that in the distributed 
burning regime, the flame width is strongly affected by the integral length scale of the 
turbulence (and hence the domain size of the present calculations). The effect of the integral 
l ength scale will b e the subject of future work, but can be predicted by Damkohler scaling 
dDamkohlerlll940h. which will be discussed in sect ion HI Some preliminary calculations (see 



Woosley et al.l (]2008al ) and lWoosley et al.l (l2008bl )) have produced encouraging results. 



4. CONCLUSIONS 

We have explored, in three dimensions, the properties of flames interacting with isotropic 
Kolmogorov turbulence for the conditions appropriate to a Type la supernova. In particular 
we have examined the consequences of a turbulent energy cascade with e = u* 3 /L* = 10 15 
erg g _1 s _1 . This might correspond, for example, to a macroscopic integral scale of 10 km 
in the supernova and speed on that scale of 100 km s . The range of densities explored is 
characteristic of the transition from the flamelet regime to distributed burning at this energy 
density. In particular, the Karlovitz number ((/l/Zg) 1 ^ 2 ) varies from 0.017 at the highest 
density to 230 at the lowest. Two regimes of burning are clearly discernible. At the highest 
density, turbulence merely wrinkles an otherwise laminar flame. The requirement that the 
flame be resolved meant that we were not able to examine a sufficiently large integral scale 
(/ ^> Ig) to see the multiply-folded flames that should be present in the flamelet regime. The 
width and speed of these individual flamelets were not greatly affected by the turbulence. 
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As the density is decreased, the flames initially remain in the flamelet regime; however, 
the increase in turbulence intensity relative to the flame speed leads to enhanced wrinkling 
of the flame. In this regime, we see significant Lewis number effects. The combination of the 
sensitivity of the nuclear reaction rates to temperature and the large Lewis number leads 
to significant variability in the local burning rate along the flame surface, with regions of 
negative curvature burning much more intensely than regions of positive curvature. 

For densities near 2.35 x 10 7 g cm" 3 (Ka pa 3), turbulence begins to tear the flame, 
altering its width and speed, but not completely disrupting the t hin region whe re burning is 



going on. This corresponds to the "thin reaction zone" regime of iPetersI (120001 ) . By 1 x 10 7 
g cm -3 (Ka = 230), however, the flame has been completely stirred and a qualitatively 
different sort of distributed burning occurs. The width of the turbulent flame is now much 
broader and it moves much faster. Turbulence takes over from radiation as the dominant 
mode of energy transport, even on small scales. Because of this, the ratio of heat diffusion 
to composition diffusion approaches unity, i.e. the Lewis number, which previously was 
very large, approaches unity. This is clearly evident in figure [111 which shows the relation 
between temperature and fuel concentration everywhere on the grid has collapsed to a line 
corresponding to advective transport only. 

This is the first time supernova flames have been simulated in the distributed regime 
in three dimensions, and our calculations may even be a first in the comb ustion commu nity 



as well. Terrestrial flames with this degree of turbulence usually go out (IPetersI 120001 ) . In 
a supernova, with its long time scale and large size, extinction is impossible until the star 
is completely disrupted. Still our results at high K arlovitz number c onfirm burning in the 



distributed regime. Assuming Damkohler scaling (iDamkohlerl 1 19401 ). the turbulent flame 



speed, St, and its width, l T , should obey the scaling relations 



l T = v^L, (6) 
for a turbulent diffusion coefficient, Dt ~ u' I, and nuclear time scale tF.. c . It follows that 



/ 7 \ 1/2 / L \ 1/2 
ST ( Ul\ ( 



sl \slIl) \t? 



nuc 



(7) 



where t^ uc is the laminar nuclear time scale. A key point is that the nuclear time scale 
is different in the laminar and distributed cases. Because of the different distributions of 
temperature and carbon abundance in the two cases (figure [TT]) . the nuclear time scale is 
almost an order of magnitude longer in the turbulent case. From the information in Table 
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1, for case (e), t^ uc = Il/sl = 6.5 x 10~ 4 s, while for the turbulent flame (figured] 
T nuc = It/ $t = 7.2 x 10~ 3 s. This gives st/sl ~ 6.4, approximately 15-20% higher than 
what was calculated (figure [7j). The overestimate is reasonable since the average turbulent 
diffusion coefficient is likely to be somewhat less than that derived for the largest possible 
length scale, i.e., Dt < u'l. 

In the distributed burning regime with Karlovitz number greater than ten and Damkohler 
number less than one, i.e. case (e), the mixing is dominated by turbulence and has a shorter 
time scale than the chemical time scale. Therefore, the overall reaction rate is limited by 
the chemistry that can occur along the burning path prescribed by figure [TTT e). and so 
th e turbulent nu c lear tim e sc ale is a constant. To gether with the calculations presented 
in I Woosley et al.l (j2008al ) and IWoosley et al.l (l2008bl ). we believe that case (e) has achieved 
the asymptotic value of t^ uc . 

For larger integral scales than those simulated here, assuming that r^ c remains fixed, 
the turbulent speed, St, and flame width, It, would both increase as as can be seen from 
equation (J7]) and the fact that v! oc I 1 / 3 . Such scaling cannot continue indefinitely though, 
without eventually encountering the limit St &u'. The flame cannot move much faster than 
the turbulence that carries it. 

For our assumed turbulent energy and composition, the length scale where Z 2 / 3 scaling 
breaks down is predicted to be 

u x 



h = {-) I, (8) 
or, for st = 1.95 x 10 4 cm s _1 at u' = 2.47 x 10 5 cm s _1 and / = 15 cm, as calculated here, 



l\ ~ 300 m. This compares reasonably well with the 155 m estimated by IWoosleyl (120071 ) 
using a very simplified representation of the flame structure. 

The actual integral scale in the supernova is much larger even than this large value. 
There, one expects a different sort of burning reflecting some distribution of turbulently 
broadened flamelets with characteristic scale l\ embedded in an overall flame brush as large 
as the actual integral scale, L* ~ 10 km. Exploratory calculations in progress suggest a 
highly variable flame width and speed in that domain that may be conducive to spontaneous 
detonation. 

As a practical consequence, the domination of turbulent transport at low density, means 
one no longer needs to resolve the laminar flame scale. It thus becomes possible to carry 
out meaningful simulations using much larger grid. In our next paper, we will demonstrate 
the validity of a simple sub-grid model for the turbulent transport and carry out three- 
dimensional calculations that take us into the unsteady burning regime. 
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Case 


(a) 


(b) 


(c) 


(d) 


(e) 


Fuel density (xl0 7 g/cm 3 ) 


8 


4 


3 


2.35 


1 


Ash density (xl0 7 g/cm 3 ) 


5.74 


2.64 


1.91 


1.43 


0.52 


Laminar flame width (Jl) 


5.00 x 10" 3 


3.19 x 10~ 2 


7.24 x 10~ 2 


1.49 x 10" 1 


2.31 


Laminar flame speed (sl) 


2.62 x 10 5 


7.52 x 10 4 


4.26 x 10 4 


2.57 x 10 4 


3.54 x 10 3 


Fuel Mach number {sl/cf) 


5.0 x 10~ 4 


1.6 x 10~ 4 


9.7 x 10~ 5 


6.1 x 10" 5 


1.0 x 10~ 5 


Laminar peak FCR 


2.10 x 10 15 


4.82 x 10 13 


8.96 x 10 12 


2.04 x 10 12 


7.61 x 10 9 


Laminar velocity change 


1.03 x 10 5 


3.85 x 10 4 


2.45 x 10 4 


1.63 x 10 4 


3.24 x 10 3 


Laminar temperature change 


4.20 x 10 9 


3.60 x 10 9 


3.37 x 10 9 


3.18 x 10 9 


2.57 x 10 9 



Table 1: Laminar flame properties. 



Case 


(a) 


(b) 


(c) 


(d) 


(e) 


Fuel density (xl0 7 g/cm 3 ) 


8 


4 


3 


2.35 


1 


Domain width (L) 


0.32 


2.0 


4.6 


9.5 


150 


Domain height (H) 


1.28 


8.0 


18.4 


38.0 


1200 


Integral length scale (/) 


0.032 


0.2 


0.46 


0.95 


15 


Turbulent intensity (u') 


3.17 x 10 4 


5.85 x 10 4 


7.72 x 10 4 


9.83 x 10 4 


2.47 x 10 5 


Gibson scale (I a) 


1.8 x 10 1 


4.3 x 10" 1 


7.7 x 10~ 2 


1.7 x 10~ 2 


4.4 x 10~ 5 


Karlovitz number ^/Il/Ig 


1.67 x 10" 2 


0.274 


0.968 


2.96 


2.28 x 10 2 



Table 2: Simulation properties. 
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Fig. 1. — Comparison of the variation of length scales at different fuel densities. The solid red 
line shows how the laminar flame width increases with decreasing density. The dashed red 
line shows how the Gibson scale (Iq) decreases with density (Kolmogorov scaling appropriate 
for a star has been assumed, see equation [Q. The flame thickness and Gibson scale are 
approximately equal at a density of p-j ~ 3. This marks the transition between the flamelet 
and distributed burning regimes. The cell width is shown by the solid blue line; the resolution 
of the turbulent simulations was chosen such that li = 4Ax. The available computational 
resources limit the size of calculation and so the domain width in each case was L = 256Ax 
and is shown by the solid black line. The forcing term used to maintain the turbulence 
imposes an integral length scale tha t is approximately a tenth of the domain width, and 
is shown by the dashed black line. I Asp den et al.l (120081 ) demonstrated that the effective 
Kolmogorov length scale using this numerical scheme is approximately 0.28 Ax, and is shown 
by the solid green line. The vertical dashed lines denote the densities of the five turbulent 
simulations. 
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Fig. 2. — Kinetic energy wav enumber spectrum f rom a homogeneous isotropic turbulence 
simulation at 256 3 , taken from lAspden et al.l (|2008l ). The same forcing term used to maintain 
the turbulence is used in the present flame calculations. This figure demonstrates the range 
of scales captured by the ILES method for turbulent flow simulations at the resolution used 
in the present study. The dashed black line denotes the minus five-thirds decay expected in 
an inertial range. 




Fig. 3. — Diagram of the simulation setup (shown in two-dimensions for clarity). The domain 
is initialized with a turbulent flow and a flame is introduced into the domain, oriented to that 
the flame propagates toward the lower boundary. The turbulence is maintained by adding a 
forcing term to the momentum equations. The top and bottom boundaries are outflow and 
solid wall boundaries, respectively. The side boundaries are periodic. 
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(a) 




-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 

Normalised height 

Fig. 4. — Normalized laminar flame profiles. Each color denotes a different quantity, as 
indicated by the legend. The line style corresponds to the five different cases; the dotted 
line is case (a) through to the solid line for case (e). The length scales have been normalized 
by the laminar flame width, defined to be the carbon width li = (AXc)/ max |VXc|, where 
Xc is the carbon mass fraction. Each quantity has been normalized by the minimum and 
maximum values attained (see table [2]), q = (q(z) — g min ) / '(<?max — <7min)- Note the velocity is 
given in the frame of reference of the flame. 
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Fig. 5. — Instantaneous vertical slices through each flame. In each case, the left-hand panel 
is fuel consumption rate, normalized by the corresponding laminar value (except for case (e), 
which was normalized by one fifth of the laminar value because it burns much less intensely 
than the laminar flame), and the right-hand panel is the temperature field, again normalized 
by the laminar value, (refer to table [2] for the values). The top legend shows the range for 
each normalized value except fuel consumption rate in case (e) , which is shown by the lower 
legend. 



Fig. 6. — Three dimensional instantaneous fuel consumption rate in each of the five cases, 
rendered such that the burning intensity is linearly related to opacity of the image. In the 
flamelet regime, the flame is coherent and flat. As the relative turbulence increases, the 
flame becomes increasingly disrupted, which leads to regions of enhanced burning and re- 
gions of local extinction. Finally, in the distributed burning regime in case (e), the character 
of the burning changes dramatically; the flame is greatly broadened, burns much less in- 
tensely (recall that the normalization is by a fifth of the laminar value), but the overall fuel 
consumption rate is enhanced. 
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Normalised time (e) 
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Normalised time (a-d) 

Fig. 7. — Normalized turbulent flame speeds. The time scale has been normalized by the 
laminar burning time, i.e. the time taken to burn one laminar flame width, specifically Il/sl- 
Cases (a)-(d) use the lower time scale; case (e), which has a significantly different evolution, 
uses the upper time scale. The turbulent flame speed has been normalized by the laminar 
flame speed. 
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Fig. 8. — Comparison of isosurfaces based on temperature (top) and carbon mass fraction 
(bottom) for case (b). The surfaces are colored by the locally integrated fuel consumption 
rate (see the text and equation [4]). The arrows labeled 'H' and 'L' highlight regions of 
high and low fuel consumption rate, respectively, both of which are correlated with negative 
curvature in the temperature isosurface. However, in the carbon isosurface, the surface is 
much lower, indicating that although the fluid is at a sufficient temperature to burn, the fuel 
is absent. 
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Fig. 9. — Joint probability density functions of the locally integrated fuel consumption 
rate (s l T ) normalised by the laminar flame speed (sl) against curvature of the temperature 
isosurface (kt) normalized by the laminar flame width (II); (a) p 7 = 8, (b) p 7 = 4, (c) p 7 = 3. 
The (negative) correlation is clear, specifically that the fuel consumption rate increases with 
negative curvature. For case (b), the regions of negative curvature labeled 'L' in figure [H 
with low fuel consumption rate, explain the bimodal PDF here. 
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Fig. 10. — Normalized turbulent flame speed per unit area. The turbulent flame speed has 
been normalized by the laminar flame speed, and the area has been taken to be the area 
of the temperature isosurface. The time scale has been normalized by the laminar burning 
time (Il/sl)- This demonstrates that a Markstein correction is required to construct a flame 
model that captures the curvature dependence of the flame speed. 
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Fig. 11. — (a-e) Joint probability density functions of carbon density against temperature. 
The red line denotes the laminar correlation. Case (a) is very close to the laminar burning 
path. Cases (b), (c) and (d) show greater variation around the laminar path. Case (e) 
present a collapse to a single curve disparate from the laminar case. This is indicative of the 
competition between turbulent mixing and thermal diffusion; in case (a) thermal diffusion 
dominates the turbulence, in case (e) turbulent mixing dominates, and in the intermediate 
cases the effects of both processes are of comparable importance. 
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Fig. 12. — First moment with respect to fuel consumption rate of the joint probability density 
function of fuel consumption rate against temperature. Again, the red line denotes the 
laminar burning path, with which case (a) is in close agreement. Case (e) has collapsed to a 
curve that is different from the laminar curve, and reiterates the lower local fuel consumption 
rate observed in this case. The intermediate cases (b) through (d) show variation around the 
mean path and move between the laminar and turbulent cases, and it is particularly clear 
that greatly increased local fuel consumption rates are observed albeit with low probability. 
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Fig. 13. — First moment with respect to fuel consumption rate of the joint probability density 
function of fuel consumption rate against carbon density for case (e); the other cases present 
similar behavior to figures tH] and [12j This figure again shows the collapse to a burning path 
different from the laminar case indicating the dominance of turbulent mixing, and the lower 
local fuel consumption rates in this case, but also shows that the burning occurs at lower 
carbon mass fractions than in the laminar case. 
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Fig. 14. — Comparison of laminar (left) and turbulent (right) profiles in case (e). In partic- 
ular, note the disparate length scales between the two figures; the turbulent flame is around 
60 times wider than the laminar flame. Each flow variable has been normalized by the min- 
imum and maximum values in the laminar flame, except for the velocity which has been 
normalized by the turbulent burning speed. Note the significant change in the shape of the 
profiles and the relative widths. In the laminar case, thermal diffusion dominates, and so 
the temperature profile (and therefore density profile) is much wider than the species mass 
fraction profiles. In the turbulent case, thermal diffusion is dominated by turbulent mixing 
and so the temperature, density and species mass fractions have similar profiles resembling 
hyperbolic tangents. Furthermore, the fuel consumption rate is greatly reduced in the tur- 
bulent case, but because of the greatly broadened flame width, the overall fuel consumption 
rate is enhanced. 



